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Abstract 

We study the effects on the dynamics of kinks due to expansions and contractions of the space. 
We show that the propagation velocity of the kink can be adiabatically tuned through slow expan- 
sions/contractions, while its width is given as a function of the velocity. We also analyze the case of 
fast expansions/contractions, where we are no longer on the adiabatic regime. In this case the kink 
moves more slowly after an expansion-contraction cycle as a consequence of loss of energy through 
radiation. All these effects are numerically studied in the nonlinear Klein-Gordon equations (both 
for the sine-Gordon and for the (p^ potential), and they are also studied within the framework 
of the collective coordinate evolution equations for the width and the center of mass of the kink. 
These collective coordinate evolution equations are obtained with a procedure that allows us to 
consider even the case of large expansions/contractions. 
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I. INTRODUCTION 



Domain walls are the interface between two regions where the system is close to two 
different stable or metastable states. Their formation and dynamics are relevant for many 
domains of Physics as Solid State Physics (magnetic materials) , Liquid Crystals Science 
P], Plasma Physics 0,0] and Cosmology (phase transitions in the Early Universe) 0, Q|- 
Spatially extended systems with more than one stable or metastable state can present domain 
walls, as is the case of systems that have undergone a first order phase transition. An 
expansion or contraction of these systems can be realized through an external action on the 
system (magnetic materials, liquid crystals, plasmas), or they are already present in the 
system dynamics (phase transitions in the Early Universe). 

The kink solutions of the one-dimensional nonlinear Klein Gordon equation for potentials 
with two or more stable states (as sine-Gordon or ip'^) provide models for the propagation of 
plane domain walls in the direction perpendicular to the domain wall. These models have 
been successfully used in several contexts as in the description of the dynamics of domain 
walls in ferrodistortive materials and the magnetohydro dynamic mode trajectory in 
reversed-field pinch experiments j^lsj . On the other hand, expansions and contractions can 
be realized in magnetic materials through the mechanical action on the material Isl, and in 
plasmas through confinement j^; while in the case of phase transitions in the Early Universe, 
the expansion is already present in the system dynamics. These expansions/contractions 
can be parameterized with a scale factor a(t) when they are homogeneous in space jioj ]. 
The characteristic times for the time variation of a{t) can be shorter or longer than the 
characteristic times of the kink dynamic. 

With these facts in mind, in this paper we consider kinks of the nonlinear Klein-Gordon 
equations for the sine-Gordon and the 0'^ potentials (as a paradigmatic examples of integrable 
and non integrable nonlinear Klein-Gordon equations, respectively), and we study theirs 
dynamics in slow and fast expanding/contracting spaces. 

In order to achieve this purpose, our paper is organized as follows: In Section lUl we 
introduce the nonlinear Klein-Gordon equations in expanding/contracting spaces, and derive 
the continuity equations for the energy and momentum densities. In the next section. Section 



im we use the Rice Ansatz 



111 ] in these continuity equations, and we reduce the problem 



with an infinite number of degrees of freedom to an approximate description in terms of two 
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collective coordinates (CC), the center and the width of the kink. In Section HV| we study 
the kink dynamic in both slow and fast expanding/contracting spaces using the collective 
coordinate evolution equations, and integrating numerically the complete nonlinear Klein- 
Gordon equations. We first consider the case of slow expansion/contraction within the 
adiabatic approximation, and later we discuss the fast expansion/contraction case and the 
departures from the adiabaticity that emerge. Finally, in the conclusions, we summarize 
and discuss the results. 

The previous sections are complemented with two appendices where some technical details 
related with expanding/ contracting spaces are reviewed, and the equivalence between several 



collective coordinates approaches [12, lid, llJ] to this problem is shown. 



II. KLEIN-GORDON EQUATION IN EXPANDING SPACES 

The Klein-Gordon equation in non expanding spaces can be obtained by extremizing the 
action defined by the Lagrangian density. 



S 



dt dx C 



dt dx 



(1) 



In an expanding (or contracting) space with one spatial dimension the space-time distance 
between events is given by 

ds'^ = dt^ - a^{t)dx^ , (2) 

where a{t) is the scale factor that gives the dilation of the physical spatial distances (see 
Appendix A). Thus, the physical spatial distance between two points with coordinates Xi 
and X2 is 

di2{t) = a{t)\x2 - xil . (3) 

The Klein-Gordon equation in this expanding space can be obtained noting that after 
the spatial dilation dx — > a{t)dx the action for a Klein-Gordon field is 



= y dtdxa{t)C = j dtdxa{t) 



-u{4>) 



(4) 



2^' 2a?{t) 

where the subscripts t and x indicate the time and space partial derivatives respectively; 
and U (0) is the nonlinear Klein-Gordon potential. 
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The Euler-Lagrange equation for the field 



d 5[a{t)C\ ^ d 5[a{t)C] 5[a{t)C] _ ^ 
dt 5(j)t dx S4>x S(f) 

gives the equation of motion, 

,„^Hm.-^^'^^0, (6) 

with H = at/a. 

We have to keep in mind that the physical distance is a{t)\x2 —xi\ and not \x2 —Xi\ (the 
latter is usually called the comoving distance). Thus, when a{t) grows the space experiences 
an elongation and when a{t) decreases the space experiences a contraction. 

We define the energy density pE{t, x), the momentum density pp{t, x), and the momentum 
current jp(t,x), as the time-time, space-time and space-space contravariant components, 
respectively, of the energy-momentum tensor. We obtain 

PE = + + (7) 
Pp = — , (8) 



JP = ^U0?+o^-W) • (9) 



1 /1^2 l-"' 

^ [-2^^ + 2 a^ 

The energy-momentum conservation (in the covariant sense) gives the following relations 

f^, ^(,,, ^.0. (10, 

^ — -t-' (U, 

that also can be obtained using the equation of motion (jUj). 

This also implies relations between the total energy E = J adx pe , the total momentum 
P = J adx pp , and the total momentum current Jp = J adx jp . They are given by 

^ + Ha'jp = 0, (12) 

[provided pp(t, oo) = pp(t, — cxd)], and 

f + 2//P^0, (13) 

[provided jp(t, oo) = j p(t, —oo)]. Note that for a space that is neither expanding nor 
contracting [i.e., a{t) = 1] E and P are conserved quantities 



III. COLLECTIVE COORDINATES ANALYSIS FOR PERTURBED KLEIN- 
GORDON EQUATION 



In the previous section we have shown that the momentum and the energy are no longer 
conserved quantities in expanding/contracting spaces, instead they satisfy the continuity 
equations ()12|) and ()13|) . respectively. We are interested in the evolution of the kink that 
we describe through two collective coordinates, its width and its center of mass position. 
Therefore, we neglect the possible energy transfer to other degrees of freedom. 

In order to obtain the equations of motion for these two collective coordinates (CCs) 
we use the Rice Ansatz in the continuity equations (fT^ and (fT^ [with the energy and 
momentum given by Eqs. ((Tj)-®], and we study two particular nonlinear potentials, the 
sine-Gordon (sG) potential U{(f)) = 1 — cos(0), and the 0^ potential U{(f)) = |(1 — 0^)^- 
The domain walls in these systems are represented by kink-like solutions that describe the 
transition between two regions in different minima of the potential. For the sine-Gordon 
potential the Rice Ansatz assumes the function 

0(x, t) = 4 arctan |exp 

whereas for the 0^, the Rice Ansatz reads 



(15) 

where X{t) and l{t) represent in both cases the center and the width of the kink, respectively, 
in the comoving frame. Note that for the unperturbed system [take in Eq. (jHl) a{t) = 1] 
the kink-like solutions are represented by Eqs. (jl4j) and (jl5|) with X{t) = vt and l{t) = 
Is = /oVl — v^, being v and Is the constant velocity and width of the kink, respectively; and 
/q = 1 (^0 = V^) the width of the soliton at rest for sG (v^^). The time dependence of the 
scale factor a{t) (expansions/contractions) implies changes in the width of the kink, as we 
show later, justifying the choice of the Rice Ansatz. 

Inserting the Ansatz, Eq. (|T^ or (|THjl . in the expressions for pp, pE, and jp and integrating 



X - X{t) 



(14) 



(x, t) = tanh 



X - X{t) 
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over X, we obtain 

= ^ — r + ^ — r + ^U W' ^ ^ 

_ Mo/o^' , aMo/o/' , Mo / /q .^^x 
•^^^^^ - ^^ + ^~^ + 2^(^^~7^J' ^^^^ 

respectively; where a = 7r^/12 and Mq = 8 for the sine-Gordon and a = (tt^ — 6)/12 and 

Mo = 2v^/3 for These expressions are replaced in the continuity equations fll2p and 

(I13p . giving the following systems of ordinary differential equations (ODE) for X{t) and l{t): 

a(t)P(t)l(t) , , 

f = -2///-. (20) 

72 l + «'Wir72 (21) 



V - 2HII - 211 

where the dots denote the derivative with respect to t. If the initial conditions are those of 



a kink with initial position X(0) and initial velocity X(0), this implies l{fi) = Iq^ 1 — ^^(O) 
and Z(0) = 0, and -P(O) is determined by Eq. (fT9|l at t = 0. The same collective coordinates 



evolution equations are obtained usin g th e Lagrangian method 



14L llfij l or the Generalized 



Travelling Wave Ansatz (GTWA) Il6l | (based on projection techniques), as t is shown 
in Appendix iBl 

This set of evolution equations, Eqs. (fTIH) - (j?H) . involves the variables X{t) and l{t) 
through the momentum P{t) defined in Section |n] Note that the equation for the mo- 
mentum is linear, and therefore it can be solved exactly yielding 



with P(0) = X(0)MoZo/ y 1 ^ X^(0). The other two equations are nonlinear and coupled, 
and they are analyzed in detail in the next section. 



IV. KINK DYNAMICS IN EXPANDING/CONTRACTING SPACES 



First of all, let us remark that the main physical variables associated with the kink 
propagation are the physical center of the kink, Xphys{t) = a{t)X{t), its physical width 
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Iphysif) = a(t)l{t), and its physical velocity 

Vphysit) = HXphys{t) + • (23) 

The first term on the r.h.s. is the contribution to the kink velocity due to the expansion 
or contraction of the space (i.e., it turns to zero when the space finishes its elongation or 
contraction H = 0). On the other hand, the second term is the peculiar velocity of the kink 
with respect to the propagating space, 

V,ecit) = ait) ^ . (24) 

For the sake of clarity, we analyze separately the cases of slow and fast expan- 
sion / contraction. 



A. Slow expansion/contraction (adiabatic regime) 

Note from Eqs. (|Tn|) . pn|l and (pT|l that the adiabatic approximation implies 

1 



a 



f - 2HII - 211 
Hence, from Eq. (j^TJ) we obtain the relation 



(25) 



lpkys{t) = lo yJl-V^Ut) 



(26) 



that links the physical width and the peculiar velocity of the kink. Furthermore, from Eqs. 
(PI) and (1221)-® we obtain 

P(0) 



Vpecit) 



(27) 



v/MoV(t) + p2(o)- 

From the previous Eqs. (j26|) and (j27|) we see that an expansion decelerates the peculiar 
motion of the kink and makes it wider, asymptotically 



lim Vpec{t) = , lim lphys{t) = k- 

a(t)— >oo a(t)— >oo 



(28) 



On the other hand, a contraction accelerates the peculiar motion of the kink and makes it 
sharper, asymptotically 



lim Vpecit) = 1 



lim Iphysit) = 0. 



(29) 
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The adiabaticity conditions can be obtained from Eq. ()25p and require only slow expan- 
sions/contractions 

m « , (30) 

pec 

where = l/{^/alo) is the so called Rice frequency for zero velocity, = vT2/7r = 1.10... 
for sG and ^r = y/6/{n'^-6) = 1.24... for 

We have studied the dynamics for different expansion and contraction rates, both for the 
sine-Gordon and the 0^ potentials. In particular. Fig. [T] (and also Fig. Ej) shows the results 
for an expansion followed by a contraction parameterized by the function 

/X . , /t-to\ Aa , /t-ti\ 

„«).l + _tanh(^^j--tanh(^^j, (32) 

with Aa, to, ti and At constants. (Aa represents the changes in a; to and ti, the times 
when the expansion and contraction take place, respectively and At is the characteristic 
time interval where the changes in a takes place.) 

In Fig.^we show the effects on a kink of a slow expansion followed by a slow contraction, 
that verify the adiabatic conditions {H ~ Aa/At = 0.1 and H ~ Aa/ (At)^ = 0.01). In this 
case, when the adiabatic conditions are verified, both the CCs equations and the adiabatic 
approximations are in good agreement with the exact results, especially for the center of the 
kink Xphysit) and its peculiar velocity Vpedf) (see Fig. H}. It is important to note that the 
agreement is good even when the variations of a are large. 

This implies that the main effects of an slow expansion/contraction are the change of 
the width and the speed of propagation of the kink, following the adiabatic relations fl2fj|l 
and (f?fj) . Therefore, these results show that the speed of a kink can he tuned by slowly 
expanding/contracting the space. 

Note that the CCs evolution equations go further and are able to predict the small 
oscillations of the kink width 



23| produced by the slow expansion/contraction (see Fig. [TJ. 
Therefore, they can compute deviations from adiabaticity excluding the radiative effects 
(because they involve the transfer of energy to other degrees of freedom). 

It is important to stress that the collective coordinate evolution equations obtained with 
an Ansatz with fixed width and variable center of the kink does not predict any variation of 
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the kink speed due to expansions/contractions. Therefore, an Ansatz that allows the width 
of the kink to evolve, like the Rice Ansatz, is an essential ingredient in order to obtain the 
correct variation of the speed of the kink under expansions/contractions. It is interesting to 
remark that, in the (h^ model the oscillations of the width are related to excitations of the 

n 

internal mode while in the sine-Gordon equation it has been shown that the excitation 
of certain phonons can imply oscillations in the shape of the kink 

On the other hand, in general the applicability of the perturbative approaches is restricted 
to slow and small expansions/contractions. The perturbative approaches rewrite Eq. (jH)) in 
the form (ptt — <Pxx + dU/dip = ef = —H{t)(f)t — {^/(^^{i) — where ef is treated as a 

small perturbation. Therefore, their range of applicability is limited not only to small H 
but also to small (1/a^ — 1), i.e., slow and small expansions/contractions. 
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B. Fast expansion/contraction (non-adiabatic regime) 

When the expansion/contraction is faster the adiabatic approximation breaks, and non- 
adiabatic effects appear as radiation (compare the kink profiles in Fig. ^ and Fig. 12)) and as 
a change in the final width and velocity of the kink [that are no longer those predicted by 
Eqs. (PBj) - ()27|) ]. See Fig. |2l that corresponds to the non-adiabatic regime H ~ Aa/At = 2 
and H ~ Aa/{At)^ = 4. The CCs evolution equations predict strong oscillations in the 
velocity and the width of the kink. In the numerical integration of the complete evolution 
equations [Eq. the oscillations are indeed present in the width. However, there are other 
degrees of freedom to which the energy can be transferred. This results in the damping of 
the oscillations through radiation emission. The loss of energy can also be easily shown 
noting that after an expansion-contraction cycle the kink moves slower (see Fig. |2)). 

We would like to stress that all the previous comments about slow and fast expan- 
sions/contractions apply both for the sine-Gordon and the 0^ potential. We have performed 
numerical simulations and comparisons of the adiabatic approximation, the CCs evolution 
equations, and the full evolution equation for both potentials obtaining analogous results. 
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FIG. 1: Sine-Gordon kink dynamics in a slow expanding and later contracting space. Top-left 
panel: scale factor a{t), Eq. with Aa = 1, to = 50, ti = 150 and At = 10. Top-right panel: 
spatial profile of the kink (f){x,tfix) at times tji^ = (solid line), tjix = 100 (long-dashed line) 
and tfix = 200 (dashed line). Three following panels: position of the center of the kink Xphys{t), 
peculiar kink velocity Vpec{t), and width of the kink lphys{t) ^ functions of time given by the 
numerical simulations of the PDE (solid line), the numerical solutions of the CCs equation (long- 
dashed line), and the adiabatic approximation (dashed line). Bottom-right panel: momentum of 
the kink P{t) exact solution. Initial conditions! 'ci kink with X(0) = 10 and X{Q) = 0.1. 
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FIG. 2: Sine-Gordon kink dynamics in a fast expanding and later contracting space. Top-left 
panel: scale factor a(t), Eq. ()32|) with Aa = 1, to = 40, ti = 100 and At = 0.5. Top-right panel: 
spatial profile of the kink (j){x,tfix) at times tji^ = 37.5 (solid line), tjix = 75 (long-dashed line) 
and tfix = 150 (dashed line). Three following panels: position of the center of the kink Xphys{t), 
peculiar kink velocity Vpec{t), and width of the kink lphys{t) ^ functions of time given by the 
numerical simulations of the PDE (solid line), the numerical solutions of the CCs equation (long- 
dashed line), and the adiabatic approximation (dashed line). Bottom-right panel: momentum of 
the kink P{t) exact solution. Initial conditions! -ci kink with X(0) = 10 and X{Q) = 0.1. 



V. CONCLUSIONS 



We have studied the effects on a l<;inl<; of slow and fast expansions/contractions of the 
media, both for the sine-Gordon and for the 0^ equations. We have used the Rice Ansatz 
jll | in the continuity equations for the momentum and the energy, in order to obtain the 
evolution equations for the width and the center of a kink in an expanding/ contracting space. 
In general this set of ODEs cannot be solved analytically and the solution must be obtained 
numerically. However, this effective description in terms of the collective coordinates usually 



gives re 



20 



21 



evant information and insight of the evolution and the phenomena involved 
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22| |. In our case the CCs approximation is able to describe very accurately the 



dynamics of the kink when the process is slow, even when the variations of a are large. We 
have found that the main effects of a slow expansion/contraction are the change of the width 
and the velocity of the kink, following the adiabatic relations fj26p and ()27|). In addition, the 
collective coordinate evolution equations for the width and the center of the kink are able to 
take into account the small deviations from adiabaticity (excluding explicitly the radiative 
effects, because they are related to the transference of energy to other degrees of freedom). 

An important consequence of the adiabatic relations is that the speed of a kink can be 
tuned by slowly expansion/contraction of space. Thus, this result provides a mechanism 
to control the speed of a kink whenever we can expand and contract the media where it 
propagates. 

On the other hand, fast expansions/contractions break the adiabatic approximation, 
giving rise to non-adiabatic effects, as for example radiation, that imply a change in the 
final width and velocity of the kink (i.e., they are no longer those predicted by the adiabatic 
relations). The collective coordinate evolution equations predict strong oscillations in the 
width of the kink. In the full problem, these oscillations are also present. However, there 
are other degrees of freedom to which the energy of these oscillations can be transferred. 
This results in emission of radiation and in the damping of these oscillations. In this sense, 
the fast expansion/contraction of space provides a method to implement a fast change in 
the speed of the kink, and also to generate phonons. 
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APPENDIX A: EXPANDING SPACES 



This appendix includes some technical details related to expanding spaces. 
1. Metric 

The metric in an homogeneous and isotropous flat space is given by 

ds^ = dt^ — a'^{t)dx'^ = g^^dx^dx^ , (Al) 
where a{t) is the scale factor for the expansion/contraction, the metric tensor is 

g,.= I ^ ° I . (A2) 
Its inverse, g^" , defined by g^vg"'^ = gj^ = 6^'^ (with a Kronecker delta) is 



/i ^ 

1 



(A3) 



The contraction with g^y lowers indices, and the contraction with g^^'^ raises indices. 

2. Invariant element of space-time volume 
Defining 

g = -Det^^, , (A4) 
it can be shown (jl^ pages 98-99) that 

y/gdtdx (A5) 

is an invariant volume element under general coordinate transformations, {dtdx transforms 
with the Jacobian, while y/g transforms with the inverse of the Jacobian.) 

Therefore the appropriate relation between the action S, and the Lagrangian density C 

is 

S = J dtdxy/gC . (A6) 

This definition implies that if the action is invariant under general coordinate transforma- 
tions the Lagrangian density is also invariant. In our case ^/g = a{t). 
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3. Energy-momentum tensor 



The energy-momentum tensor is 



Some useful relations to calculate its components are 



(A7) 



9> = g''^d^(j) = g^^d. 



H ■ 



(A8) 
(A9) 



The time-time component of the energy-momentum tensor or energy density is 

k2 



PE = r« = - g' 
1,0 



2a2 



the space-time component or momentum density is 

Pp = T^* = T*^ = 9*05^ 
and the space-space component or momentum current is 



(AlO) 



gt^C = -Mx , 



(All) 



(A12) 



APPENDIX B: EQUIVALENCE AMONG DIFFERENT CC APPROACHES 

In this appendix we show that the evolution equations for the collective coordinates 
obtained with the derivation of the momentum and the energy (that we used in Section 



and with the so- 



III) are the same as those obtained with the Lagrangian method 

called CTWA HQ (based o„ p.o.eetion teeh„i,ues). In all eases, we use the R,ee Ansa. 
as an approximated kink-like solution of the following perturbed nonlinear Klein-Gordon 
equation 



a2(t) 



(Bl) 



where Pit) is a time dependent damping coefficient and z{x,t,(f)) represents a generic per- 
turbation on the system. Note that for [3 = and z{x,t, 0) = we recover the system ®, 
introduced in Section II. 
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1. Generalized Lagrangian Formalism 



Introducing a new time variable 

dt' 



c(t 



(B2) 



where c(t) is given by 



c{t) = expj^ [H{t') + (3{t')]dt' 
a{t) 



«(0) uo 
Eq. ()B1|) becomes the dissipationless equation, 

(t>TT (J)XX dU 



exp{ I (3{t')dt'} , (B3) 



C2(r) A2(r) 

with C(r) = c(t(r)), v4(r) = a(t(r)), and Z{x, r, (p) = z{x, t{r), (j) 
The evolution equation ()B4|) can be obtained from the action, 



+ Z(x,r,0), (B4) 



S = j dTdxC = J dTdxS^^<l>l-^^^(l)l-C^iT)[Ui<l>) + Mix,T,<l>,<l),)]'j (B5) 

C being the Lagrangian density, and M{x, r, (p, (j)^) an "Euler-Lagrange integral" of 
r, 0). Some of the possible integrations are 

• the (p functional integral, 

M{x,T,(f)) = - j Z{x,T,(p) , (B6) 

• and the following integration for the particular case where Z{x, r, 0) = F{x, t)G{4>) 

M(x, r, 0, 0,) = 0, G(0) f dx'F{x',T). (B7) 



XO 



Once we define the Lagrangian, 

L = I dxC , (B8) 



the Euler-Lagrange evolution equations for the collective variables are obtained inserting a 
given Ansatz in this Lagrangian and using either the expression (|B6|) or ()B7|) . This procedure 
is the generalization of the Legrangian formalism developed in jl^ 

15 



Using the Rice Ansatz defined by ()14|) and psp (for the sG and (p^ potentials, respectively) 
and after some straightforward calculations, we obtain the Lagrangian as a function of our 
collectives variables and theirs derivatives with respect to r 

L(X,X',M') = ^{Xf + ^ilf- (B9) 
Replacing L{X, X' ,1,1') into the Euler-Lagrange equations, 

(BIO) 
(Bll) 



d 1 




dL 




M'J 


~ dX 


d 




dL 


dr 


[dFJ 





and rewriting the equations of motion in the time variable t, we get 



. a{t)Pmt) 

Mo/o ' ^^''^ 
dP _ „ 1 

'di 



-{(3 + 2H)P-^y^ dez{X + ei,t,(l))(l)eY (B13) 



a 



21 



2 



/ + 00 
dez{x + ei,t,(f))e(f)e, (bi4) 
'OO 



Mo/o 

where 6 = [x — X{t)]/l and the dots denote the derivative with respect to t. These are 
the evolution equations for the collective variables obtained for either of the two previous 
expressions of M{x,t,(f),(j)x) E<1- f)B6|) or Eq. ()B7|) (in their respective regime of validity). 
Taking in these equations /5 = and z = 0, we obtain the Eqs. (fTI?|) . (pUj) and (j?T|) derived 
in Section III. 



2. GTWA 

In order to apply the GTWA we rewrite Eq. (IB1|) as 

(t) = i), 

a^[t) d(j) 

The procedure to obtain the CCs equations with the GTWA consists of inserting our specific 
functional form for 0, Eq. ()14j) for sine-Gordon or (fT3j) for 0^ potential, into Eq. ()B15|) . mul- 
tiplying the first equation by dip/dX and the second one by d(j)/dX, taking their difference 
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and integrating over x, equating the result with zero; and repeating the same procedure with 
dip/dl and d(j)/dl. This gives to the same ordinary differential equations for and 
l{t) as those obtained in the previous subsection [see Eqs. ()B12|) - fjB14|) ]. 
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